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ABSTRACT 

We examine the angular momentum transport properties of disks composed of 
macroscopic particles whose velocity dispersions are externally enhanced ("stirred"). 
Our simple Boltzmann equation model serves as an analogy for unmagnetized fluid 
disks in which turbulence may be driven by thermal convection. We show that inter- 
particle collisions in particle disks play the same role as fluctuating pressure forces and 
viscous dissipation in turbulent disks: both transfer energy in random motions associ- 
ated with one direction to those associated with another, and convert kinetic energy 
into heat. The direction of angular momentum transport in stirred particle and fluid 
disks is determined by the direction of external stirring and by the properties of the 
collision term in the Boltzmann equation (or its analogue in the fluid problem). In 
particular, our model problem yields inward transport for vertically or radially stirred 
disks, provided collisions are suitably inelastic; the transport is outwards in the elastic 
limit. Numerical simulations of hydrodynamic turbulence driven by thermal convection 
find inward transport; this requires that fluctuating pressure forces do little to no work, 
and is analogous to an externally stirred particle disk in which collisions are highly 
inelastic. 

1. Introduction 

When turbulence in unmagnetized Keplerian disks is driven by thermal convection, angular 
momentum appears to be transported inwards (Ryu & Goodman 1992; Kley, Papaloizou, &; Lin 
1993; Stone & Balbus 1996; Cabot 1996). If generically true, this implies that convection cannot 
provide the "enhanced" angular momentum transport (relative to collisional viscosity) needed to 
account for observations of accreting systems in astrophysics. 

Despite appearances, inward transport does not contradict the fact that an isolated Keplerian 
disk relaxes to its minimum energy state by transporting angular momentum outwards (Lynden-Bell 
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k, Pringle 1974). The reason is that simulated convective disks are not completely isolated. They 
tap a source of free energy other than their intrinsic differential rotation: an externally imposed, 
superadiabatic entropy gradient. 

A more mechanistic understanding of the behavior of turbulent fluid disks can be obtained by 
considering disks composed of discrete macroscopic particles. Transport properties and stability 
criteria obtained under one picture often carry over to the other, to within dimensionless factors of 
order unity. The fluid/particle dualism has been employed with considerable success in studies of 
planetary rings (see, e.g., Goldreich & Tremaine 1982) and galactic structure (see, e.g., Chapter 6 
of Binney & Tremaine 1987). 

Indeed, one-to-one correspondences can be established between terms in the Reynolds stress 
tensor equation describing turbulent fluid disks, and terms in the Boltzmann stress tensor equation 
describing particle disks. This paper explores the transport properties of the latter system, with the 
aim of shedding light on the former. We demonstrate how simple particle disk models can reproduce 
many of the transport phenomena found in numerical simulations of vertically or radially convective 
fluid disks. The particle viewpoint also provides insight into the stability properties of isolated, 
unmagnetized, fluid disks. An analysis similar to ours was carried out by Kumar, Narayan, & Loeb 
(1995); this is discussed in §5. 

In §2, we write down the Boltzmann stress tensor equation governing particle disks and dis- 
cuss its intimate relationship with the Reynolds stress tensor equation governing hydrodynamic 
turbulence in a fluid disk. In §3, we prove a general theorem that relates the direction of angular 
momentum transport to the behavior of the velocity dispersion tensor in the presence of external 
driving. Section 4 is the most important of this paper; there we consider a simple model for the 
collision term in the Boltzmann equation that allows one to easily solve for the stress tensor of 
the disk and, we believe, capture much of the relevant physics regarding both stability and the 
direction of transport. In §5, we summarize and discuss our results. 



The second moments of the Boltzmann equation describe the evolution of the viscous stress 
tensor. For an axisymmetric disk, in cylindrical coordinates, these read (see, e.g., Borderies, Gol- 
dreich, & Tremaine 1983): 



2. Stress Tensor Equations for Particle Disks 



dtPrr ~ ^Pr<f> = (%>rr) c + H rr 
K 2 

dtP4>4> + -^Pr<t> = {dtPu) c + H H> 
dtPzz = (dtPzz) c + H zz 



(1) 



(2) 



(3) 




(4) 
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where we have neglected third moments. In equations ([|)-@, ^ is the angular frequency, k 2 = 
r~ 3 d(Qr 2 ) 2 / dr is the square of the epicyclic frequency, and pij are the components of the stress 
tensor. The terms (dtPij) c describe changes in velocity dispersions due to interparticle collisions. 
Collisions perform three tasks: they (1) extract energy from the mean shear into random motions, 
(2) transfer energy in random motions associated with one direction to those associated with an- 
other, and (3) convert kinetic energy into heat. The terms Hy are inserted to simulate external 
stirring of the disk, with iJ r< £ = so as not to explicitly induce angular momentum transport. 



2.1. Boltzmann vs. Reynolds 

Since the fluid equations are derived from the Boltzmann equation it is not surprising that 
hydrodynamic turbulence in fluid disks satisfies energy conservation equations very similar to equa- 
tions (|l])-(||). The Reynolds stress equation can be obtained by taking the Navier-Stokes equation, 
decomposing the velocity and pressure into mean and fluctuating components, and multiplying the 
entire equation by the velocity. This yields (e.g., Speziale 1991; Kato & Yoshizawa 1997): 

dtTij = —TikdkUj — TjkdkUi + ir^ — €ij — dkCijk, (5) 

where U is the mean velocity of the fluid, u is the fluctuating velocity, Tij = (u{Uj) is the stress 
tensor, and ( ) denotes an ensemble average. The quantity is a sum of turbulent and viscous 
fluxes, 

Kij = -(mdjp + Ujdip) (6) 

and 

eij = 2v(d k Uid k Uj), (7) 

where p is the fluctuating gas pressure and v is the kinematic viscosity. The first two terms on the 
right hand side of equation (|5j) describe how turbulence extracts energy from the shear in the mean 
flow. The terms nij are known as the "pressure-strain" correlation in the fluid dynamics literature;^ 
the diagonal terms of TTij are simply the work done by fluctuating pressure forces associated with 
the turbulence; among their other tasks, these terms attempt to enforce isotropy of the turbulence 
by transferring energy from velocity fluctuations associated with one direction to those associated 
with another. The terms represent the viscous dissipation of energy on small scales (note that 
en is positive definite). 

It is straightforward to evaluate the above expression for dtTij in cylindrical coordinates. We 



4 To be precise, our iiij is related to the usual pressure-strain correlation tensor by a total divergence, which we 
have absorbed into CV,fc. 
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neglect the flux divergence^] to obtain 

K 2 

dtTtj,,)) + -rr T r( p = — £00 (9) 

d t T zz = tt zz - e zz (10) 

K 2 

dtT r( p + 2Q^rr — 2J7T00 = 7r r — e r< ^ . (11) 

Equations (|8[) — (p~0|) are identical to the energy equations given by, e.g., Balbus & Hawley (1998; see 
their equations 58-60). 

Equations (0)-@ for the evolution of the stress tensor in a particle disk are very similar to 
equations (P)-(|IT|) for the evolution of the stress tensor in a fluid disk. In particular, the rotational 
profile of the disk enters identically in the two cases; both f2 and k 2 couple to the stress tensor 
components in exactly the same way for particle and fluid disks. Moreover, the "source" and "sink" 
terms in the energy equations are analogous; the collision term, (dtPij) , in the Boltzmann equation 
plays the same role as — 6™ in the Reynolds stress equations. 

This comparison indicates that the dynamical role of epicycles is the same in both fluid and 
particle disks, and that any difference in the angular momentum transport properties of the two 
disks can only be due to differences in the sources and sinks in the energy equation. We demon- 
strate this formally in the next section by proving a theorem that relates the direction of angular 
momentum transport to the properties of the collisional and stirring terms. This theorem is very 
useful for interpreting the model problem of §4. 



3. The "Long Axis Must Get Longer" Theorem 

Theorem: For thin axisymmetric disks in which d£l/dr < 0, steady inward transport of angular 
momentum requires that the net tendency of collisions and external stirring be to further lengthen 
the principal equatorial long axis of the velocity ellipsoid and to further shorten the principal 
equatorial short axis. This theorem applies to fluid disks as well as to particle disks so long as the 
interparticle collision term, (dtPij) c , is replaced by TTy — eij. 

Proof: We transform equations (|l])-(||) to the principal axis basis of the velocity ellipsoid. Let 
pu and P22 be the greater and lesser of the two components of the stress tensor, respectively, in 
the equatorial plane: 

Pll = [Prr +PW+ \J (Prr ~ Pcjxj,) 2 + 4p^ ] / 2 (12) 
P22 = [Prr + PW ~ \J (Prr ~ PwY + 4 Pr<t> 1 / 2 > ( 13 ) 



5 This assumes "homogeneous" turbulence, i.e., that the flux varies on a length scale ~ r, much greater than the 
vertical scale height of the disk. 
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where we shall always take the positive square root. By symmetry about the midplane, p 3 3 = p zz . 
The orientation of the velocity ellipsoid is given by 5, the angle between the longer equatorial axis 
and the radial direction: 



tan 5 = — 



P11 - 

2p r <t> 



(Prr -P<t>4>) + \J {Prr ~ Pcj>c/>) 2 + ^P 



2 

r<j> 



(14) 



valid for \5\ < ir/2 and nonisotropic Pij^\ Since the denominator in equation (14) is positive, 

sgnp r(f) = sgn<5. (15) 

Outward (inward) transport corresponds to positive (negative) p r< f„ which in turn corresponds to 
positive (negative) 5. 

In the principal axis basis, the viscous stress equations ([l])-(|]) read (cf. Goldreich & Tremaine 
1978): 

d tPll - Sn (sin 25) p u = (d tPll ) c + H n (16) 
dtP22 + Sn (sin 25) P22 = (d tP22 ) c + H22 (17) 

%>33 = {dtP33) c + H 33 (18) 

(pn-p 22 )d t 5-(2-S S m 2 5)np 22 + (2-Sco S 2 5)n Pll = (p n - p 22 ) (d t 5) c + H 12 , (19) 

where S = — d\nn/d\xir. 

Equations (^?|) and ( p!7[ ) state that for S > 0, the steady (dt = 0) inward transport (5 < 0) of 
angular momentum demands that the net tendency of collisions and external stirring be to further 
lengthen the principal long axis of the velocity ellipsoid and to further shorten the principal short 
axis. 

Collisions between hard, imperfectly elastic spheres tend to isotropize the velocity ellipsoid: 
(dtPu)c < < (dtp 22 ) c - Thus, in the absence of stirring, transport is normally outwards in 
axisymmetric particle disks. 



4. Angular Momentum Transport in Externally Stirred Disks 

In this section we solve for the steady state structure of a particle disk using a simple model 
for the effects of interparticle collisions and external stirring. We also discuss the fluid analogue of 
our particle disk model in some detail. 



If S — ±7r/2 or pij is isotropic, p r( j> = and there is no transport. 
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4.1. The Particle Scattering Function 

Following, e.g., Cook & Franklin (1964), Shu & Stewart (1985), and Narayan, Loeb, & Ku- 
mar (1994) we take the scattering function in the Boltzmann equation to be given by the Krook 
approximation : 

(dtf)c = (20) 

T 

where r is the collisional mean free time, / is the particle distribution function, and fo is the 
"equilibrium" or "post-scattering" distribution function. The corresponding collision term in the 
viscous stress equations is (dtPij) c = (p®j —Pij)/i~. The physical idea behind the Krook model is 
that when particles collide they are removed from phase space at a rate —f/r and are returned to 
phase space with the new distribution function _fo- 

We assume that collisions do not generate off diagonal stresses so that = o if i ^ j. We 
account for the fact that particle collisions are inelastic, converting some fraction of the kinetic 
energy of the colliding particles into heat. With this model, the relevant second moments of /o are 
given by (see Kumar et al. 1995) 

2 2 2 

P° rr = (1 " £r) y , P% = (1 - , and p° zz = (1 - ^ , (21) 

where a is the velocity dispersion (a 2 = p rr + p^ + p Z z) and the parameters £j e [0, 1] measure the 
inelasticity of the collisions; £ ~ is elastic, while £ ~ 1 is highly inelastic. 

Equation (^Tj) allows the inelasticity of colliding particles to be anisotropic, i.e., to be different 
for the r, eft, and z directions. Such an anisotropic scattering function is unrealistic for real particle 
disks in nature (e.g., planetary rings); collisions between hard spheres at the microscopic level 
cannot be expected to respect the cylindrical coordinate system. Our strategy is to first solve for 
the structure of "real" externally stirred particle disks, in which the £j's are equal. This captures 
much of the physics in which we are interested. We will show, however, that the assumption of 
equal £j's does not furnish a fully adequate model of fluid disks and that the introduction of an 
"anisotropic inelasticity" can remedy this defect. 



4.2. The Fluid Analogue of the Scattering Function 

Our particle scattering function has a mathematically rigorous analogue in fluid turbulence. 
Namely, it is equivalent to a particular closure scheme for the a priori unknown turbulent quantities 
€ij and iTij: 

e i0 = ^T and Try = ^^ 2 <%- (22) 

Interpreted as such, r is roughly the timescale on which free turbulence would decay (the eddy 
turnover time) and the inelasticity parameter, £j, measures the efficiency with which fluctuating 
pressure forces perform work in the i—th direction. The closure model proposed in equation ( p2[ ) 
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differs from conventional second-order closure models of the Reynolds stress equation (e.g., Speziale 
1991; see Kato & Yoshizawa 1997 for an application to Keplerian disks). We will show, however, that 
this simple model provides considerable insight into the angular momentum transport properties 
of accretion disks. 

In contrast to r, it is not obvious what values to assign £j in hydrodynamic turbulence. Al- 
though in one sense turbulence is "inelastic," with energy cascading to small scales in the absence of 
source terms, the physical meaning of defined by equation (|22|) is more precise: fluid turbulence 
is "inelastic" if and only if fluctuating pressure forces are inefficient at doing work. Throughout 
the remainder of this paper it is worth keeping in mind that numerical simulations suggest that 
the highly "inelastic" limit is appropriate for Keplerian disks. Note also that, unlike in the particle 
disk case, we would not expect the &'s to be identical in the fluid case. The fluctuating pressure 
forces are macroscopic in origin and thus can plausibly be different for the r, <j>, and z directions. 



4.3. The Stress Tensor Equations 

With our model for the interparticle collision term, the steady state equations for the stress 
tensor become 

p rr - (1 - £ r )y - 40rp r< ^ = H rr , (23) 
2 2 

+ ^Pr<p = H H , (24) 

2 

Pzz-(l-£*)Y=flzz, (25) 



and 

~^Prr ~ 20rp^^ + p r(t> = 0, (26) 
where H^j = rHij, i.e., the net stirring over the mean free time. 

An "energy equation" can be derived by summing equations (|23|)-(p5|): 



H = H rr + H H + H zz = io 1 + 2 Pr<t> n , T, (27) 



where Q' = dft/dlnr and 3£ = £ r + £^ + £ z . For 17' < 0, equation (27) shows that, if there is no 
external stirring, transport must be outwards or zero in steady state. This is because any energy 
lost due to inelasticity must be extracted from the shear in the mean flow, which requires p r( f, > 0. 
No such thermodynamic constraint on the direction of angular momentum transport exists if there 
is external stirring. 
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4.4. Transport in a Keplerian Disk 



The general solution to equations (|23|)-(|26|) is given in the Appendix. Here we focus on the 
case of a Keplerian disk, for which k = 0; in the next subsection we briefly contrast the Keplerian 
disk with a constant angular momentum disk. We define x = Sir and solve for the r(f> component 
of the stress tensor: 



Pr<f> = 



xa 



/ j [1 - £d - I [1 - - H rr j + AH H { \ 

\ 1 + x 2 [2.hH rr + 10^0 + 4# 22 ] J 



(28) 



where Ha = Ha/H. The velocity dispersion of the disk is not an independent parameter but is 
determined by solving the equation 



a 



H + x 2 (2.5H rr + 10H H + 4H ZZ ) 
(1 + Ax 2 )^ - 0.5£ r x 2 + 2^x 2 - 1.5x 2 ' 



(29) 



4-4- 1- Real Particle Disks: Isotropic Inelasticity 

For a real particle disk with isotropic inelasticity (£r = = £z = £) equations ( p8[ ) and ( p9| 
become 

xa 2 ( i - e - H„i + 4# 00 e 



1 + x 2 



2.5H rr + 10H M + 4H Z 



(30) 



and 

^2 = H + x 2 (2.5H rr + lOff 00 + 4H ZZ ) 

a (1 + 5.5x 2 )£ - 1.5x 2 ' { ' 

In the limit of no external stirring, we recover the usual results of particle disk theory. In such 
an isolated disk, the inelasticity adjusts to the optical depth of the disk, 



«-TT^- (32> 

1 + 5. oar 



This follows from equation (|3l]) by setting H = and requiring finite a 2 . Since the collisional 
inelasticity depends on the relative impact velocity, equation (|32|) determines a velocity dispersion 
(Goldreich & Tremaine 1978); call this a®. Inserting equation (|32] ) for £ in equation (|30|), we confirm 
that the transport is always outwards, with[| 

£(J 2 xo~ 2 

^ = ^ = 2TTT^- (33) 



7 To derive eq. write eq. [g(J as p r< ^ = 0.5xa 2 N/D. Rewritting the denominator as D — 1 + ix 2 + x 2 (6H$^ 
l.5H rr ) and substituting eq. |32] for f into the numerator, N, all terms oc Hij cancel, yielding eq. fall. 
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The fact that isolated fluid disks appear to be hydrodynamically stable can be explained 
mechanistically in the sense that they do not satisfy equation (32). For any x, equation (pJ3) has a 



solution only if £ e [0,0.27]. Simulations suggest that the work done by fluctuating pressure forces 
is far less efficient than this. Therefore, for H = 0, the only solution to equation (p_T|) is ctq = 0, 
i.e., no turbulence. 

What happens if the disk is externally stirred? For H <C £oq the disk structure is only weakly 
affected by the stirring; corrections to p r( f,, e.g., are 0[H/£(Tq] <C 1. Physically, this is because, 
for o"o 7^ 0, heating must become significant relative to the "natural" velocity dispersion in order 
to modify the disk structure. For a hydrodynamically stable fluid disk with do = 0, however, 
the "natural" velocity dispersion vanishes; thus, for any H, the properties of the external stirring 
(together with the "collisions" ) determine the angular momentum transport properties of the disk. 

In the presence of significant stirring, i.e., H <; Oq, the direction of angular momentum trans- 
port can reverse; in particular, the sign of the angular momentum transport is then determined by 
the sign of the numerator in equation (pjl). Note that the sign of p r( f, does not explicitly depend 
on x = rO. This is consistent with the arguments given in §2 and §3; only the properties of the 
collisions and the stirring determine the direction of angular momentum transport, not the degree 
to which the disk is collisional. 



From the numerator in equation (|30| ) it follows that (1) energy input along r tends to lead 
to inward transport (H rr enters with a negative sign), (2) energy input along tends to lead to 
outward transport (H^ enters with a positive sign), and (3) elastic scattering tends to lead to 
outward transport (p r< j, > if £ <C 1). All of these results are readily understood using the theorem 
of §3. The principal long axis of the velocity ellipsoid generally has a much larger projection along 
the radial direction than the azimuthal direction.^] Thus r stirring corresponds to elongating the 
long axis of the velocity ellipsoid and so should (and does) lead to inward transport (and vice- 
versa for (p stirring). The transport only respects the direction of external stirring, however, if 
the collisions are fairly inelastic. In the elastic limit, collisions efficiently transfer energy from one 
direction to another, isotropizing the distribution function; transport is then always outwards. 

Two special cases of equation (|3(]) are of particular interest: 

1. R stirring only: This is relevant for interpreting the work of Stone, Pringle, & Begelman 

(1999) and Igumenshchev & Abramowicz (1999) who carried out simulations of advection- 
dominated accretion flows, which are radially convective. In these flows, one can treat the 
convection as "externally" driven by the viscous dissipation associated with magnetic fields. 
As discussed by Narayan, Igumenshchev, & Abramowicz (2000) and Quataert & Gruzinov 

(2000) , the convective angular momentum transport is inwards in these simulations, with 
a as large as pa —0.1 (a pa Pr^/cr 2 is the dimensionless viscosity). As noted above, it is 



8 In the collisionless unstirred limit, equation ( |l9[ ) states that p r r/p<j><p = 2/(2 — S) = 4 for steady-state Keplerian 
disks. 
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particularly natural for radial convection to produce inward transport since radial convection 
tends to directly expand the long axis of the velocity ellipsoid. More concretely, equation ( |30"| ) 
predicts that the transport should be inwards in the presence of r stirring so long as £ > 1/2. 
Moreover, for £ ~ 1 and x ~ lj^| a should be ~ —0.1 in the presence of r stirring, consistent 
with simulations. 

2. Z stirring only: This is relevant for interpreting the simulations of Stone & Balbus (1996) in 
which vertical convection in an otherwise stable Keplerian disk generated low levels of inward 



angular momentum transport (a ~ — 10~ 4 ). For a particle disk, equation (30) shows that 
transport is always outwards in the presence of z stirring: p r< ^ oc (1 — £). This can be under- 
stood by noting that the z dynamics in thin reflection-symmetric disks is largely decoupled 
from the r and (j) dynamics which determine the direction of transport. In particular, energy 
input in the z direction is transferred to the r and eft directions by particle collisions, which 
attempt to isotropize the distribution function. By the theorem of §3 the transport must then 
be outwards. 



4-4-%- "Anisotropic Inelasticity" : Inward Transport with Vertical Stirring 

The analysis of the previous section applies to real particle disks in which the inelasticity 
is isotropic. Interestingly, it shows that vertically stirred particle disks would transport angular 
momentum outwards, opposite to what is found in simulations of hydrodynamic turbulence. As 
noted in §4.2, the fluid analogue of inelasticity is inefficient fluctuating pressure forces; such forces 
are macroscopic in origin and therefore need not be the same in the r, <f>, and z directions. We now 
discuss how such "anisotropic inelasticity" can produce inward transport in the presence of vertical 
stirring. 

In the presence of vertical stirring only, equation ( |28| ) reduces to 

xa 2 f j[l-^]-|[l-£r] 
Vr4, ~ 2 [ l + 4x 2 

From the numerator in equation ( |34| ) it follows that elasticity in the r direction tends to lead to 
inward transport (1 — £ r has a negative prefactor), while elasticity in the (j) direction tends to lead 
to outward transport (1 — £^ has a positive prefactor). 

Inward transport in the presence of z stirring requires > 0.75 + 0.25£ r . Physically, this 
means that fluctuating pressure forces are particularly inefficient at doing work in the eft direction 
(cf. Stone & Balbus 1996). In fact, in the absence of explicit 4> stirring, equation ( |28| ) shows that 




9 This means that the eddy turnover time is comparable to the rotational period of the disk, as is indeed the case. 
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the transport is always inwards if £^ — ► l.[^| In order to reproduce a ~ -10" 4 in the presence of 
z stirring, as found in the simulations, equation ( |34[ ) requires (taking x ~ 1) 1 — £ r ~ 10~ 3 , and 
1 - ^ < 10" 3 . 



4.5. Transport in a Constant Angular Momentum Disk 

Balbus and collaborators have emphasized the interesting point that a constant angular mo- 
mentum, k = 0, disk is analogous to a linear shear flow (e.g., Balbus, Hawley, & Stone 1996). 
Indeed, numerical simulations show that the k = disk is unstable to nonlinear perturbations 
and produces outwards angular momentum transport. This feature is also readily captured by our 
model. By equation (|A1| ) the r<f> stress in a k = disk is given by 

which is positive definite. For unstirred k = disks in which the particle inelasticity is isotropic, 

8x 2 



«=3T^' < 36 > 



In contrast to equation (32) for a Keplerian disk, equation ( p6| ) can be satisfied for some x for any 
£ e [0,1]. For a fluid disk this means that, regardless of the efficiency of fluctuating pressure forces, 
outward transport can be sustained in an unstirred disk. This is consistent with the fact that shear 
flows are nonlinearly unstable. 



5. Discussion 

This paper explores the angular momentum transport properties of axisymmetric, externally 
stirred particle disks as a means of better understanding hydrodynamic turbulence in fluid disks. We 
are not, however, advocating that turbulence is a sort of macroscopic kinetic theory, with colliding 
turbulent "blobs" literally replacing colliding particles. Rather, we argue that it is physically 
revealing to exploit the similar mathematics governing turbulent fluid disks and particle disks. 

There are one-to-one correspondences between terms in the Boltzmann stress tensor equation 
governing particle disks and terms in the Reynolds stress tensor equation governing fluid disks (§2). 
This analogy reveals that the role of differential rotation and epicyclic motion is the same in both 
fluid and particle disks; both the rotation rate and the epicyclic frequency couple to the components 
of the stress tensor identically in the two cases. Consequently, the transport properties of fluid disks 



1 This is analogous to Stone & Balbus's argument that axisymmetric simulations of hydrodynamic turbulence will 
always find inward transport. 
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can differ from those of particle disks only to the extent that the "scattering function" for fluid 
elements differs from that of particles. The formal, and very useful, analogy is that interparticle 
collisions in particle disks play the same role as viscous dissipation and fluctuating pressure forces 
in turbulent disks. Both transfer energy in random motions associated with one direction to those 
associated with another, and convert kinetic energy into heat. 

Motivated by the above considerations, we solve for the structure of a particle disk using a 
simple model for the collision term in the Boltzmann equation (§4). We allow different collisional 
inelasticities along the r, and z directions. Although an anisotropic inelasticity is not realistic for 
actual particle collisions, its introduction serves to make the analogy with fluid turbulence complete. 
The fluid analogue of particle inelasticity is the inefficacy with which fluctuating pressure forces 
do work, and there is no a priori reason for the latter to be isotropic. A final ingredient of our 
model is the insertion of an external stirring mechanism, modeled as constant source terms in the 
r, <j>, and z energy equations. External stirring is relevant for interpreting numerical simulations in 
which either vertical or radial convection is present (e.g., Stone & Balbus 1996; Stone, Pringle, & 
Begelman 1999). 

We believe that this particle model captures the transport properties of fluid disks. With a 
single equation for the stress tensor (eq. []28|| ) we can reproduce results from standard particle disk 
theory and from a number of numerical simulations of hydrodynamic turbulence. 

In the limit of no external stirring, differential rotation is the only source of free energy in 
our disks. The angular momentum transport induced by turbulence or particle collisions must 
then be outwards or zero in steady state. This is because both particle collisions and turbulence 
generate heat, which can only be compensated for by a positive rcf> stress (see eq. ^7|). To have 
outward transport in a Keplerian disk requires that "collisions" be relatively efficient at isotropizing 



the distribution function (see eq. [32| and the associated discussion). This tendency towards 
isotropization is normally satisfied in axisymmetric particle disks such as planetary rings; these 
systems therefore transport momentum outwards while maintaining some degree of random motion 
(turbulence). By contrast, the hydrodynamic stability of isolated fluid disks can be understood 
by noting that, according to simulations, fluctuating pressure forces are extremely inefficient at 
isotropization. Thus there can be no self-sustaining turbulence in these systems.^ 

In the presence of significant external stirring, there are no thermodynamic constraints on the 
direction of angular momentum transport. Instead, the sign of the transport is determined by 
the nature of the stirring and the properties of the collision term in the Boltzmann equation (or 
its suitably defined analogue in the fluid problem); formally, a necessary and sufficient condition 
for inward transport is that the combined tendency of collisions and stirring be to make the long 



n Kato & Yoshizawa (1997) have argued that if pressure fluctuations in turbulent Keplerian disks behave like 
they do in shear flows and other turbulent laboratory fluids, the isotropization is sufficient to yield self-sustaining 
turbulence. This assumes, however, that Reynold's stress closure models developed for laboratory problems are 
applicable to Keplerian disks. 
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equatorial axis of the velocity ellipsoid yet longer and the short equatorial axis yet shorter (§3). 

More physically, for any mechanism of external stirring (e.g., convection, spiral shocks, etc.), 
our model predicts the direction of transport as a function of the "inelasticity" of collisions (eq. p8[ ). 
The latter is, of course, not known a priori for the fluid problem, which underscores the parametrized 
approach of our model. Nonetheless, it is gratifying that, so long as collisions are suitably inelastic, 
our model finds inward transport for both vertical and radial stirring (e.g., convection), with the 
transport naturally much more efficient for the latter than for the former. Inward transport in 
the presence of purely vertical stirring requires, however, an anistropic "inelasticity," for which 
"collisions" are more inelastic in the azimuthal direction than in the radial direction (§4.4.2). By 
equating inelasticity with inefficient fluctuating pressure forces, as advocated above, we find results 
that are consistent with analytic arguments and simulations of hydrodynamic turbulence (e.g., 
Balbus & Hawley 1998). 

To conclude, we note that our analysis expands on several aspects of the Boltzmann equation 
model for radial convection in rotating media considered by Kumar, Narayan, & Loeb (1995). 
These authors treated "collisions" between convective eddies using the same scattering function 
we have considered (eq. pi~| ) an d found that for sufficiently inelastic collisions, radial convection 
would produce inward angular momentum transport. Interestingly, we now see that this inward 
transport has little to do with the details of convectively unstable fluids (e.g., the superadiabatic 
entropy gradient, correlated temperature and velocity fluctuations, etc. analyzed by Kumar et al.). 
Rather, it is a generic property of disks in the presence of external "stirring" of the radial velocity 
dispersion. Moreover, we have emphasized the analogy between inelastic collisions in particle disks 
and inefficient fluctuating pressure forces in turbulent disks, better establishing the relevance of 
kinetic theory models to the fluid problem. 

We thank Steve Balbus, Roger Blandford, Peter Goldreich, Jeremy Goodman, Pawan Kumar, 
Ramesh Narayan, Scott Tremaine, and an anonymous referee for useful discussions and comments. 
EQ is supported by NASA through Chandra Fellowship PF9-10008, awarded by the Chandra X-ray 
Center, which is operated by the Smithsonian Astrophysical Observatory for NASA under contract 
NAS 8-39073. EC acknowledges partial support from a Caltech Kingsley Foundation Fellowship. 
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A. Solution for General k 



Denning x = fir, y = kt, and x' = dx/dlnr = (y 2 — Ax 2 ) /2x, one can solve equations (|23| 
to find: 



Prtf, = xa 



§[i-fc]-&[i-e 



H K 



and 



a 2 = 



1 + 4y 2 + ^H rr - Axx'H, 



H + 4Hy 2 + £fH rr - 4xx'H 4 



(l + VX-fg(l-^)-4,(l-^) 



(Al) 



(A2) 
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